Asymptotic dynamics of coined quantum walks on percolation graphs 
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Quantum walks obey unitary dynamics: they form closed quantum systems. The system becomes 
open if the walk suffers from imperfections represented as missing links on the underlying basic 
graph structure, described by dynamical percolation. Openness of the system's dynamics creates 
decoherence, leading to strong mixing. We present a method to analytically solve the asymptotic 
dynamics of coined, percolated quantum walks for a general graph structure. For the case of a circle 
and a linear graph we derive the explicit form of the asymptotic states. We find that a rich variety 
of asymptotic evolutions occur: not only the fully mixed state, but other stationary states, stable 
periodic and quasi-periodic oscillations can emerge, depending on the coin operator, the initial state 
and the topology of the underlying graph. 

PACS numbers: 03.67.Ac, 05.40.Fb 



Quantum walks (QWs) are natural generalizations of 
random walks, obeying unitary quantum evolution 
The possible applications of quantum walks range from 
quantum information to modeling transport phenomena 
0. Quantum walks Q, similarly to classical random 
walks, can serve as a tool to design efficient algorithms 
At the same time, their conceptually simple defini- 
tion with only slight modifications makes them a versatile 
model for a variety of physical phenomena, motivating in- 
tensive research to understand their properties in recent 
years. In a series of very recent experiments Q several 
aspects of quantum walks have been demonstrated. The 
next generation of quantum walk experiments aims at 
demonstrating possible fields of application of walks in 
solid state physics @ or open quantum dynamics Q. 

Randomly removing edges from a graph one arrives at 
the problem of percolation [8j . Percolation is perhaps the 
simplest, but a very generic model for disordered media. 
In static percolation, the set of randomly missing edges 
is fixed throughout the propagation, whereas in dynamic 
percolation some random noise can change the connec- 
tions among the vertices in each time-step @. In the 
context of quantum walks, randomly broken links are a 
natural source of noise: in the quantum information set- 
ting they would correspond to errors, whereas in trans- 
port models they are a consequence of random variations 
of the medium. 

The effect of broken links in coined quantum walks was 
first studied in 10] as a general source of decoherence. 
The long time behavior of the diffusion coefficient for the 
percolated walk on the line was calculated in [ll|. and its 
time evolution was numerically investigated in [12| . The 
asymptotic limit for the global chirality distribution of 
a quantum walk on the line with percolation restricted 
to the half-line was considered in [131 ] . Static percolation 
for coined quantum walks was studied in the context of 
spreading 12] and search efficiency 2M an d, furthermore, 
for continuous-time quantum walks in [l5| . 

Time evolution of discrete time quantum walks is uni- 



tary by definition, it can be specified by the coin and the 
step operators. The iteratively defined evolution is then 
realized by repeated application of the coin and the step 
operator maintaining deterministic evolution of the den- 
sity operator. Randomness can enter the quantum walk 
either by applying measurement [Ril - fl8j or introducing 
an external random variable. External disturbance may 
affect only the coin or only the step operator. Fluctua- 
tions in the coin operator usually lead to decoherence and 
the transition from ballistic to diffusive motion. While a 
lot of interest was paid to changes of the coin [19( much 
less is known about varying the step operator. Percola- 
tion can be viewed as a perturbation of the step operator. 
Other types of perturbation also lead to interesting direc- 
tions, e.g. the realization of Levy flights in QWs where 
connections with a distant vertex can app ear, fundamen- 
tally modifying spreading properties 20] . 

In this Letter we resolve the asymptotic limits for 
coined quantum walks on finite graphs with dynamical 
percolation. We obtain explicit analytic forms for the 
asymptotic evolution for walks on a finite linear graph 
and prove that generally there exist initial states from 
which the completely mixed state cannot be reached. 
We construct the attractor space of the superoperator 
governing the evolution of the density operator. Beside 
various steady states, we find stable asymptotic periodic 
oscillations as well as quasi-periodic behavior. 

Let us consider d-regular simple graphs G(V,E) with 
a finite number N = \V\ of vertices. We introduce dy- 
namical bond percolation on this family of graphs as fol- 
lows: Every edge in this graph has the same independent 
p probability of being present, and the complementary 
1—p probability of being broken (missing) . In every step 
of the quantum walk we randomly choose a configura- 
tion of the edges K, C E. On each of these configurations 
we define a unitary quantum walk by combining the step 
operator and the coin operator. We assume the coin op- 
erator to be always the same. Hence the definition of the 
step operator is alternated. At positions with missing 
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edges we invoke the concept of reflection. 

The Hilbert-space of discrete time coined quantum 
walks (QWs) on percolated graphs defined above is the 
tensor product of the Hp position space and the He coin 
space. Hp is & ./V-dimensional Hilbert-space spanned by 
orthonormal basis vectors corresponding to the vertices 
of the graph, and He is d-dimensional, spanned by coin 
basis states where d is the regularity of the graph. The 
coin state indicates the direction of the nearest neighbor 
vertex the walker is displaced during the next time evo- 
lution step. Let us use the notation \a) ® \b) = \a, b), 
where \a) G Hp and \b) e He- 

The unitary time evolution of the QW on a given per- 
colated graph is defined as 

U K = S K ■ (Ip <8 C) , (1) 

where 

Stc = E E \a®d,d){a,d\ 

a,d \(a,a®d)i£K 

+ (Ip®R)\a,d){a,d\ \ , (2) 

(a,a®d)<£K J 

and C is the so called coin operator from the U (d) group. 
By a © d we denote the next neighbor of vertex a in the 
direction d, and R (reflection) is in general a traceless 
permutation matrix ensuring the unitarity of Sic- In the 
examples considered in the present paper, R further sim- 
plifies to a simple transposition matrix. The effect of S/c 
could be considered as conditional displacement if the ac- 
tual edge is present and reflection if the edge is broken. 

The randomness is introduced in the process by choos- 
ing different edge configurations K. at each step, and thus 
random application of different unitary step operators. 
We denote the probability of choosing a configuration K., 
and corresponding unitary steps U/c by ttk.{p)- Since this 
process involves classical randomness, the state of the 
walker is best described by a density operator p(n) , after 
n steps. The complete state of the system changes during 
time evolution according to 

p(n + 1) = MpWWUt = * (p(n)) , (3) 

K 

where we introduced the linear superoperator 

Numerical calculation of the superoperator in Eq. ([3]) 
over all possible \K\ = y / 2 d ' N edge configurations by 
brute force methods would have an exponential compu- 
tational cost as a function of the number of vertices N. 
In the following we show that this cost can be reduced, 
requiring resources scaling polynomially with N. From 
now on we use the short notation for the matrix ele- 
ments Xf'j = (s, c\X\t, d). Utilizing Eq. © for matrix 
elements, and changing the sum over elements with the 



sum over configurations we arrive at the expression 

#S(*+D= E &)\Y,^{pW>iu^A . 

a,b.q.r \K.CE J 

(4) 

The second sum is over all possible configurations. The 
number of configurations \K\ grows exponentially with 
N . Studying matrix elements of U/c — S/c ■ C, we rec- 
ognize that only the elements of neighboring vertices are 
nonzero: Sic is sparse. Consequently, the seemingly ex- 
pensive second sum could be taken just on edges between 
neighboring sets of vertices £, thus JC C £ C E. Since 
every vertex has d neighbors, all Sjc operators map a 
single vertex to one of the d neighboring vertices. Thus 
a single run of the second sum, restricted to the set £, 
contains only 2 2d = 4 d additions in the worst case. The 
first summation is 0(N 2 ) (polynomially) dependent on 
the number of vertices, thus computation of a single time 
evolution step is reduced to the polynomial regime with 
respect to N. Moreover, the superoperator $ is also a 
sparse operator, and the cost of its construction is re- 
duced to the polynomial regime as well. 

Decoherence caused by percolation leads to mixing of 
states, and one might naively expect that asymptotically 
it reaches the completely mixed state on a finite graph. 
Examining a particular system, however, reveals that this 
is not the case. The joint position-coin limit distribution 
is not necessarily uniform, one can obtain (quasi) periodic 
behavior, as shown in Fig. [TJ 

The asymptotic dynamics of the map $ given by © 
can be described via spectral decomposition [2l|. The 
process is somewhat similar to finding the eigenvalues 
and eigenvectors of a matrix, but physically the resulting 
eigenoperators are not necessarily valid density opera- 
tors. First, we have to find all possible solutions of the 
equations 

U K X x ,iUj- = XX \ i , with |A| = 1 , for all K C E , (5) 

such that they form a maximal orthonormal set {X\^}. 
The subspace spanned by {X\,i} is called the attractor 
space, which includes density operators. The index i dis- 
tinguishes different mutually orthonormal solutions with 
a certain eigenvalue A. The density operators, contained 
in the attractor space, are the possible asymptotic states 
of the system. The next task is to take into account the 
initial condition, represented by the density operator po 
at time t = 0. In terms of the attractors, the asymptotic 
evolution (n> 1) is then given by 

/}(„) = $«G5o)= E ^{xl.p Q )X kti . (6) 

i,|A|=l 

In case the attractor space contains more than one valid 
density matrices representing different physical states, 
then asymptotic cycles may occur. Depending on the 
size and structure of the attractor space, the asymp- 
totic evolution either converges to a single state or leads 
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FIG. 1. Distance of the joint coin-position distributions from 
the uniform distribution A4(p, r) = Yl \Pi ~ r i\ f° r percolated 
quantum walks in dependence on the number of iterations. 
Initial states are parametrized according to Eq. (|16|1 , a) 
Hadamard walk on the 7-cycle with 6 — tt/2, (j> = —tt/2. Inset: 
quantum fidelity with respect to the totally mixed state, b) 
The same walk as in a) on the 8-cycle. The steady state here is 
not totally mixed, c) 7-linear graph with 9 = —tt/2, <j) = — tt. 
An asymptotic limit cycle occurs, d) Special walk on the 
8-cycle using the coin C(n/2, y/2) (see Eq. (|11|) ') with the 
same initial state as in c) . The asymptotic behavior becomes 
quasi-periodic. In all cases p = 0.5, the points are connected 
for visual aid. 



to a periodic or a quasi-periodic limit cycle. As it was 
proven in (2ll | , the generalized eigenvectors corresponding 
to eigenvalues A with |A| = 1 are all in fact eigenvectors 
(the corresponding Jordan blocks are one-dimensional). 
Moreover, the phases of these eigenvalues determine the 
eventual (quasi) periodicity of the corresponding asymp- 
totic cycle. Examining Eq. (JSJ) , one can see that it con- 
tains an equation for every unitary operator Uk.. These 
operators, in turn, are determined by the possible graphs 
with nonzero probability. Except the extremal cases of 
p = or p = 1, all graph configurations occur with finite 
probability, thus the set of unitaries is always the same 
for p greater than and smaller than 1, therefore the 
asymptotic dynamics is independent of p for < p < 1 . 

Employing ([TJ we rewrite the set of eigenvalue equa- 
tions ([5J) into 



(Ip®C) X x (I P ®C r ) 



(7) 



which allows us to find the attractors in two steps. Eq. 
([7]) represents a set of equations where the left hand side 
is always the same. First we have to find a subspace 
of matrices solving conditions imposed by different shift 
operators S^X^Sic = S^ k ,X\Sk.' for all pairs of different 
edge configurations K,,K! . Once we restrict the set of 
attractors to the latter subspace, transitivity implies that 
it is sufficient to resolve an equation of type (0 only for 
a single chosen edge configuration. The most favorable 



choice is the graph where all edges are missing, resulting 
in the condition 



(I P ® RC) X X (I P ® C^R 1 ") = XX x 



(8) 



Using a block form X\ = J2 s . \s)(t\ eg) X^'^ we rewrite 
© into the set of A*" 2 block equations of the same form 
(RC) X (s ^ (RC) ] = XX is ^\ which can be turned into 
an eigenvalue problem 



(RC) ® (RC)*x (s <^ = Xx^ 



(9) 



The vector x^'^ is defined via the relation (a, 6|aj( s,t ') = 
(a|X( s ' 4 ) |6). The importance of equation (|9|) is twofold. 
First, it shows that the spectrum of the operator RC 
determines all possible candidates of the attractor spec- 
trum. Second, eigenvectors corresponding to a given 
eigenvalue A of the operator (RC) <g> (RC)* reveal con- 
ditions imposed by the coin operator C on each block 
matrix X^ s '^ of a possible attractor X\. 

The just described method in principle allows for solv- 
ing the percolation problem for an arbitrary graph. The 
form of the solution indicates the possibility to have non- 
stationary asymptotics (the eigenvalues can have arbi- 
trary phases while being of unit magnitude). This is a 
particularly interesting result in view of the general ten- 
dency of open systems to have a unique, single asymp- 
totic state. To illustrate the power of our method we give 
the analysis of the simplest nontrivial 1 dimensional case. 

We apply the proposed method to solve the dynamical 
percolation problem on two topologically different finite 
graphs: the linear graph and the circle, both consisting 
of N vertices. The linear graph can be viewed as a circle 
with a single constantly broken link thus realizing reflect- 
ing boundary conditions, or vice versa the circle can be 
viewed as a line with periodic boundary conditions. The 
step operator takes the form 

& = EE E \aQ(~l) d ,d)(a,d\ 



d=0,l \(a,aQ(-l) d )GK 



E 

(a,ae(-l) d )<ZK 



\a, 1 — d) (a, d\ 



(10) 



Here the operations © and denote modulo-A^ oper- 
ations in case of the cycle and ordinary operations plus 
and minus in case of the linear graph. This reflects differ- 
ent boundary conditions which take place in definitions 
of both graphs. The reflection operator is chosen R = a x . 
Let us consider a two parameter class of coins 



ie la sin (3 cos /3 
cos f3 ie ia sin (3 



(11) 



with f3 tt/2. The case f3 = tt/2 leeds to degeneracy 
[22| . This class of coin operators includes the Hadamard 
coin for a = tt/2,/3 — 7r/4. Examining the structure of 
Eq. ([5]), one can first determine the spectrum of the 
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FIG. 2. Position mixing time (circles) for a percolated 
Hadamard walk on the 7-cycle as a function of distance 
threshold e. Initially, the walker was localized on a vertex 
with totally mixed coin state. The line represents our ana- 
lytic estimation Eq. QlSp . 



matrix a x C(a, j3), which has two different eigenvalues 
e ±I/3 . Consequently, the possible solutions of Eq. ©, 
i.e. attractor spectrum, can only be a subset of the set 
{l,e ±2 ^}. 

Following the recipe outlined after Eq. ([7]), after some 
tedious but straightforward calculation one can explicitly 
determine the basis of the attractor space. For the linear 
graph (reflecting boundary conditions) the dimension of 
the attractor space corresponding to eigenvalue A = e 2l/3 
is one, its base reads 



X 



t,d 



1 

2N 



exp (iaS)(— 1) 



t+d 



(12) 



For the conjugate eigenvalue, A = e~ 2l/3 , the base can be 
chosen as Y — . The attractor space corresponding to 
the third eigenvalue, A = 1, is three dimensional, its base 
{Z 1 , Z 2 , Z 3 } can be chosen in the form 

Z 1 =I 2N /V2N, (13) 
(Z 2 Y'l = exp [ia(S - 1)](1 - (-1) S )/2V2N , (14) 



>t,d 

73\ S - C 

J )t,d 



exp (ia5)(l + (- 



(15) 



with S = s — t — c + d. 

In the case of the cycle (periodic boundary conditions) 
we have additional constrains and the dimension of the 
attractor space will depend on the length of the graph. 
If N = 2k + 1 and a = ln/N {I € IN) the attrac- 
tor space has dimension equal to 2, all attractors cor- 
respond to eigenvalue A = 1 and its base can be chosen 
as {Z\1/V2(Z 2 + e la ( N ~^Z 3 )}. When N = 2k and 
a = 21tt/N, the structure of the attractor space is the 
same as for the percolated linear graph, for N = 2 and 
a =/= In the attractor space is two-dimensional, all attrac- 
tors correspond to A = 1 and its orthonormal base can be 
chosen as {Z 1 ^ 3 }. For all other cases on the cycle we 
found that there is only one trivial solution of Z 1 . These 
results clearly specify the topology induced difference be- 
tween the walk on the percolated linear graph and circle. 
This method becomes far more involved when going to 
higher dimensions. However, also in this case the method 
offers analytic results [22|]. Let us illustrate the results 
obtained for coin C(a, 0). We classify the initial states of 
the walk starting from one given position |xo), with re- 
spect to the asymptotic state it reaches on a linear graph 



or on cycles with 5 dimensional attractor space. We use 
the standard Bloch-sphere notation 

h = \xo){xo\®(^\$oMM + \( 1 - p ) I c \ > ( 16 ) 

where \4>e,<t>) = cos||0) + e 1 ^ sin §|1). Examining the at- 
tractor space, it is straightforward to see that the only 
initial coin state leading to a single stationary asymp- 
totic state, which is not completely mixed, corresponds 
to 9 = tt/2 and P > 0. The asymptotic steady state 
will be completely mixed only if the initial coin state is 
totally mixed (P = 0). The purity of the asymptotic 
state (Tr (p 2 as ) = £|A|=i \Tr{X x p } | 2 ) can be explicitly 
calculated for localized initial states, yielding 



Tr(/5 2 s 



1 

2N 



P 2 
2N 2 



(17) 



We immediately see that the most pure asymptotic state 
can be reached from a pure (P = 1) initial coin state. 

The onset of the asymptotic regime can be character- 
ized by the mixing time [23| . One can employ the fol- 
lowing definition for the position mixing time t m i X (e) — 
ml{t:M \p P (t'),p P (oo)] < e, V*' > t}, where p P (t) is 
the position probability distribution of a QW at time t, 
and AA(p, r) = £ |p, — rt \ is the Manhattan distance. If 
one assumes that the slowest decaying states correspond 
to the largest |A| ^ 1 eigenvalue of the superoperator $ 
(in the simplest case there is only one such eigenvalue), 
then one can estimate the mixing time 



loge/log|A'|-log|0|/log|A'| 



(18) 



where A' is the highest magnitude |A| ^ 1 eigenvalue of 
superoperator <I>. By O we denote the overlap of the 
subspace corresponding to A' with the initial state of 

the QW: = Tr (x\,p(0)\. We illustrate this result 
in Fig. H 

In summary, we presented an explicit method to solve 
the asymptotic dynamics of percolated QWs. Unlike in 
typical environment induced decoherence processes, in 
the considered case the asymptotic open system dynam- 
ics can in general be oscillatory. We emphasize that the 
asymptotic evolution does not depend on the percolation 
rate p if < p < 1. Already the analysis of the simple 
case of a linear-type graph provides nontrivial conclu- 
sions. First, the possibilities for the asymptotic evolution 
on the linear graph are fundamentally different from that 
of the cycle. This is a consequence of the topology of the 
graph. One can expect the occurrence of similar, but 
possibly more involved topological effects for more com- 
plex graphs, e.g. for higher dimensional lattices. Second, 
periodic boundary conditions lead to sensitivity to the 
number of sites. 

In all the discussed cases, the position distribution 
was asymptotically uniform, and asymptotic dynamics 
was limited to the chirality degree of freedom. One can 
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imagine that more general graphs open the way for non- 
uniform asymptotic position distributions. 
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